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Abstract 

Let fj.be a, finite positive Borel measure on [—1, 1], m a fixed natural number and £( Q ' /3 ) [/] — 
(1 — x 2 )f" + (/3 — a~ (a+f3+2)x)f' , with a, f3 > —1. We study algebraic and analytic properties 
of the sequence of monic polynomials (Q„)„> m that satisfy the orthogonality relations 

j £ (a ' l3) [Q n ]{x)x k d[i{x) = for all < k < n - 1. 

A fluid dynamics model for source points location of a flow of an incompressible fluid with 
preassigned stagnation points is also considered. 
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differential operators, hydrodynamic 
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1 Introduction 

Let P be the space of all polynomials. We say that a sequence {p n }n£i, + C P, is a polynomial 
system if for each n the nth polynomial p n is of degree n. 

Let \x be a finite positive Borel measure on [—1, 1] and {i n }nez + the corresponding system of 
monic orthogonal polynomials; i.e. L n (z) = x n + . . . and 

(L n ,L k )^ = J L n (x)L k {x)dfj,{x) | ^ p ifri^fc' ( L1 ) 
Denote by £( Q >^) the Jacobi differential operator on the space P, with a, (3 > —1, where 

C (a,0) [f] = il _ x 2 )f „ + (p_ a _( a + p + 2)x)f , t f e p 5 (L2) 

or equivalents (cf. [TU1 (4.2.2)]) 

L [/J " (1-^(1 + ^ ' /eF ' 
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From (|1.2I) it follows that / and C^ a '^ [/] are polynomials of the same degree. It is straightfor- 
ward that integrating (|1.3I) with respect to the (a, /3)-Jacobi measure dp a ^{x) = (1 — x) a (l+x)' 3 dx 
on [—1, 1], we obtain 

[ [f](x)dt, a j(x)=0, /GP. (1.4) 

We say that Q n is the nth monic orthogonal polynomial with respect to the pair (£( a:/3 ',/i) if 
deg[Q„] < n and 

J £ ( - a '^[Q n ](x)x k dn(x) =0 for all < Jfe < n - 1. (1.5) 

This nonstandard orthogonality with respect to a differential operator was introduced in [1], 
where the authors prove conditions of existence and uniqueness. The starting points of this work 
are [HE], where the orthogonality with respect to the differential operator C$ [/] = / + (z — £)/', 
( £ C, was analyzed. 

From (jl.ip . we have that a monic polynomial Q„ of degree n is orthogonal with respect to 
(C {a ^\p) if and only if it is a polynomial solution of the differential equation 

C^[Q n ]=X n L n , where \ n = = -n(l + n + a + p). (1.6) 

As will be shown, it is not always possible to guarantee the existence of a system of polynomials 
{Qn}neZ+ orthogonal with respect to the pair (C^ a '^ , fx). Let to G N be fixed, a fundamental 
role in this paper is played by the class V m (a,/3) defined as the finite positive Borel measures p, 
supported on [—1,1] for which there exist a non negative polynomial p of degree to, such that 
dp(x) = p^ 1 (x)dp a ^(x). 

This manuscript deals with some algebraic and analytic aspects of the sequence of orthogonal 
polynomials {Q n }n>m orthogonal with respect to the pair {L^ a, ^\p). We provide asymptotic 
results for the sequence of the orthogonal polynomials with respect to (£( a 'P\ p) and study the set 
of accumulation points of their zeros as well. In particular, we prove 

Theorem 1.1. Let p G V m {c*,/3), where to G N. If {Cn}n>m is a sequence of complex numbers 
with limit £ G C \ [— 1, 1] and {Q n } n >m the sequence of monic orthogonal polynomials with respect 
to the pair (C^ a '^\p) such that Q n {Cn) = 0, then the accumulation points of zeros of {Q n } n>m 
are located on the set E = £ (Q) 1J[— 1, 1], where £{Q) is the ellipse 

£{Q :={z£C:z = cosh(r? c + i6),0 < 9 < 2tt} , (1.7) 
andr, c := In |v?(C)i = In |C + VC 2 ~ 1|- If HO = inf |C - %\ > 2 then E = £ (£). 

— 1<X<1 

Let <p(z) — z + \J z 2 — 1 be the function which maps the complement of [—1, 1] onto the exterior 
oftheunit circle, where we take the branch of V z 2 — 1 for which \<p(z)\ > 1 whenever z G C\[— 1,1]. 
If p G V m (<X)l3), let vt,i/2, • • • j Vm G C be the to zeros of the polynomial p(z) — r Ilti( z ~ v i)i f° r 
wich dp{x) = p~ 1 (x)dp at p(x). 

For all z G C \ [—1, 1] we define the function $(p, z) and a constant <f> m as 

^ n?) - ¥>(^) \27r7 Vl -t 2 J 

The function (f> m z) is a particular case of the well known Szego's function (cf. [B] §6.1]). 
For the sequence of monic orthogonal polynomials with respect to the pair [C^ a ^\p) the 
following asymptotic behavior holds 
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Theorem 1.2. Let {Cn}n>m be a sequence of complex numbers with limit ( e C \ [—1, 1], m G N, 
/i G "P m (a,/3) ; and (Q n ) n >m the sequence of monic orthogonal polynomials with respect to the pair 
(C (a ^\fi) such that Q n (Cn) = 0, then: 

1. Uniformly on compact subsets of = {z G C : |<^(.z)| > I'/KOI}; (i- e - the exterior of the 
ellipse E(C) ) 

(U '- ] —r <t> 2 m *(p,z). (1.8) 



Pn ( z ) n^oo 

2. Uniformly on compact subsets of = {z G C : |</>(V)| < 1^(01} \ [ — 1> 1 

Qn(z) 



(i.9) 



If 5(C) > 2 i/ien (|1.9|) /io/<is /or f2 = {z G C : |y(^)| < 1^(01} (*- e - the interior of the ellipse 

£«))■ 

The paper is organized as follows. In Section[2]we study the existence of a system of polynomials 
{Qn}n>m, orthogonal with respect to the pair {C^ a ^\^). Section [3] is devoted to the study of 
recurrence relations and location of zeros of the polynomials {Q' n } n >m- In Sections |4] and [5] we 
study the asymptotic behavior of the polynomials {Q„}„> m , {Q n }n>m and its zeros respectively. 
In the last section, we introduce a fluid dynamics model for the interpretation of the critical points 
of Q n - 

2 Existence and uniqueness 

It is known that the operator has a system of monic eigenpolynomials {i» } n ez + and a 

sequence of constant (A„)„ e z + (eigenvalues), such that 

£(°^[p( a '«]=A n F^«, (2.1) 

where the eigenpolynomial Pn"'^ is the nth monic orthogonal polynomial with respect to the 
(a, /3)-measure of Jacobi d(x a ,p(x) — (1 — x ) a (1 + X Y dx on [—1,1], i.e. 

(p^\x k ^ ^ = J P^ a ' e \x)x k dp, a ^{x) = Q, for all < fc < n- 1. (2.2) 

Furthermore, from [ID] (4.21.6) and (4.3.3)] 



p( a ,0) p(a,/8)\ 

± n i ± n n 

T(n + a + 1) T(n + f3 + 1) T{n + 1) T(n + a + (3 + 1) 
2 _( 2n+ a+^+i) p( 2n + a + P + 2) T(2n + a + (3 + 1) ' 



(2.3) 



Theorem 2.1. Let n be a fixed natural number and [i a finite positive Borel measure on [—1,1]. 
Then, the differential equation (|1.6p has a monic polynomial solution Q n of degree n, which is 
unique up to an additive constant, if and only if 



L n {x)dn a ,fs{x) = 0, (2.4) 



where L n is as (|1.1[) . 
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Proof. 

Suppose that there exists a polynomial Q n of degree n, such that £( Q ' /3 ) [Q„] = — n(l+n + a + ft) L n , 
where L„ is the nth monic orthogonal polynomial for u. Then (|2.4I) is straightforward form (|1.4p . 

Conversely, suppose that L n , the nth monic orthogonal polynomial with respect to /i satisfies 
(|2.4|) . Let Q„ be the polynomial of degree n defined by 



n-l 



? n (z) = P^.«(z) + £ a„, fc P fc (a '«(*), (2.5) 

fe=0 



where a„ ; o is an arbitrary constant and 

a„,fc = v^&n.fc = t-^- ( L n,Pj: a '^) > fc = l,...,n-l. 



From the linearity of C^ a '^ [•] and (J2HJ) we get that C^'^ [Q n ] = -n(l + n + a + ft) L n . □ 

We are interested in discussing systems of polynomials such that for all n sufficiently large they 
are solutions of ()1.6[) . In this sense, the next corollary is fundamental. 



Corollary 2.2. Let n be a finite positive Borel measure supported on [—1, 1] such that d/i(x) — 
r(x)d/j, a ,p(x), where r 6 L 2 (fi a ^). Then, m is the smallest natural number such that for each 
n > m there exists a monic polynomial Q n of degree n, unique up to an additive constant and 
orthogonal with respect to (C^ a, ^\ fi), if and only if r~ x is a polynomial of degree m. 

Proof. Suppose that m is the smallest natural number such that for each n > m there exists a 
monic polynomial Q n of degree n, unique up to an additive constant and orthogonal with respect 
to (£( Q ' /3 \/i). According to Theorem [Q 



-^—jL n (x)d[i(x) = J L„ (x)dfj, a ^(x) tn = Z ' 

But this is equivalent to saying that — — = > CkL^ix) with c m ^ 0. The converse is straightfor- 

r(x) 
v ; fc=0 

ward. □ 

From the previous corollary, if fi S V m {a, ft) then the differential equation (jl.6p has an unique 
monic polynomial solution Q n of degree n for all n > m, except for an additive constant. 

Let {Cn}n>m be a sequence of complex numbers, where m € N is fixed, and assume that 
/i G V m (a,ft). We complement the definition of the sequence {Q n }n>m in (|1.5|) by considering 
that henceforth Q n for each n > m is the polynomial solution of the initial value problem 

£(.<*,&) [ y ] = \ n L n , n>m, r „„s 

y(Cn) = o. {2 ' b) 

We say that {Q n }n>mi is the sequence of monic orthogonal polynomials with respect to the 
pair (£( a 'P\ n) such that Q„(C„) = 0. 

If Q n is the monic polynomial of degree n defined by the formula 



Tn—k 



k=0 

then, the initial value problem (|2.6|) has the unique polynomial solution 

= Q n (z) = Q n (z) - Q n (C„). (2.8) 
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3 The polynomial Q' n 

Let m £ N, {CnjneN a sequence of complex numbers, and /i € V m (ct, j5) be fixed, then for all n > to 
the polynomials Q„ (solution of (|2.6p ) are uniquely determined by (|2.7p ~ (|2.8p . Without loss of 
generality, wc will complete the sequence of polynomials Q n for all n € N as follows 



Qn(z) = {pi a ' p Kz)-pi a ' p Hc n )) 

min(m.n) 

+A„ £ ( P »^ )(z) - ,» > 1 (3-1) 

k—1 

Qo(z) = l 

For convenience, only in the previous formula we consider Ao = 1. Note that {Q n }neN defined by 
(13. ip is a system of polynomials, such that Q n ((n) = for all n > 1. Let us remark that if n < to, 
in general, £ (a ^ ) [Q„] ^ A„i„. 

Additionally, as the degree of a polynomial is invariant under the operator £( a '^[-] and the 
polynomial Q n , for all n < m, is of degree n (see (I3.ip ). 

(1, [Qi] , . . . , [Q m ] , L m+X , . . . , L„, . . .) (3.2) 

is a polynomial system. 

In the following result, we show that for n > (2to + 1) the derivatives of the system of polyno- 
mials Q n satisfy a recurrence relation with a fixed finite number of terms. 

Theorem 3.1. Let \x 6 T > m (a,/3) where m G N. Then if R is any primitive of p, for each 
n > (2m + 1) the sequence of polynomials Q' n satisfies the relation 

m+l 

R(z)Q' n (z)= Y 0R,n,n-kQ'n-k( z )> ( 3 - 3 ) 
k— — rn— 1 

where the initial values Q' m+ i, ■ ■ ■ , Q'2 m +2 are given by the derivatives of (|2.7I) and 

9R,n,n-k = T (A n eR. n ,n-k + dn,n-fe)> 

eR,n,n-k = 7^— / R{x)L n - k {x)L n {x)dn(x), (3.4) 
k = f Ll(x)dn{x), (3.5) 
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hi-k . ... 
3=3i{k) 

ji(fc) = max{- and j 2(h) = min{m + 1, m + k}. 

J4(fc) 



Cn,n—k 



hn,n—j Gn—j,j — k 



J=33(k) 



33$) = max{0, k — 1} and ji(k) = min{m, k + 1}, 
c n ,i = —n, 

2n (a - (3) (n + a + (3 + 1) 



Cn.O 



(2n + a + p) (2n + a + /3 + 2) ' 

(n + a) (n + /3) (n + a + 0) (n + a + [3 + 1) 



(2n + a + (3) 2 ((2n + a + f3f - 1) 
Before starting the proof of Theorem 13. 11 we will state and prove some lemmas. 
Lemma 3.2. Let p e V m (a, ft) where m G N. ITien /or n> m 



L n (z) = ][>,„-* P£f(*), (3.6) 



/c=0 



p(z)p(-^)(z) = r „£^L„ +fc( ,), (3.7) 
where bij = ^- (^Li, P^ a ^^ , Tj as in (|2.3[) and ij as in (|3.5[) . 

Proof. As /i € V m (a, (3) then (L n ,x k ) a = for all n > m and k < n — m. Hence, o n! fc = for 
fe = 0, 1, . . . , n — to — 1 and (|3.6[) is established. 

From the orthogonality relations and (|2.2I) . if i < j or i > j + to, we have 

L i (x)Pj a ' 0) (x)dp a ^(x) = J L i (x)P^ ) (x)p(x)dp(x) = 0. (3.8) 

The relation (|3.7[) is straightforward from the Fourier expansion of pPn*'^' in terms of the (Lfc), 
fc = 0.1 n + m and (I3T51). □ 



From the lemma above, the polynomials Q n defined in (|2.7|) and its derivatives can be written 
as a linear combination of the polynomials {£ ra }nez + as we show in the next lema. 

Lemma 3.3. Under the conditions of Lemma \3.2l for n > to the polynomials Q n satisfy the 
following relations: 

Q n (x)x k dp a} p(x) = 0, fc = 0, 1, . . . , n — m — 1, (3-9) 

m+1 



(1 - z 2 )p(z)Q' n (z) = d "-k,k Ln-k(z). (3.10) 

Proof. From (|2~7| the relations (|3l)]) and 



fc=— m— 1 



6 



(i - z 2 ) Q ' n (z) = (i - z 2 ) g;(.) = a„ £ (! - - 2 ) ( p i-f (*)) . (3-n) 

fe=0 " k 

follows directly. 

Using the structure relation fulfilled by Jacobi polynomials (see QUI (4.5.5)-(4.5.6)]), we have 

(1 " ^ = Cn-Kl PttiM + Cn- k ,0 P^[Z) + C M 

Substituting this formula into ([3. lip , we obtain 

m+l 

(l-Z 2 )Q' n (z)= 5,,n-^-fW> 



k=-l 



and from (I3.7[) . (|3.10p immediately follows. □ 

Proof of Theorem \3.1[ As the sequence {Q n }nez+ is a system of polynomials, then the sequence 
of its derivatives {Q' n }neN is also system of polynomials. Hence, the polynomial RQ' n can be 
expanded as linear combination of the polynomials {Q' n }ntEN, i.e. there exist (n + m) constants 

&R,n,l, • ■ • , @R,n,n+m SUCh that 

n-1 

R(z)Q' n (z) = e R>n ^ k Q' n _ k {z). (3.12) 

k— — m 

Let £( a ' /3 ) be the linear differential operator on the space of all polynomials P defined by 
£(«,«[/']=£(«,/?)[/] for all/ GP, i.e. 

£{a,fl [/] = (i _ x * )f > + _ a _ (Q + £ + 2 )x)f. 

Since {L n }nez + is a system of polynomials and £^ a,/3 -'[-] is a linear operator, the polynomial 
[i?Q£J can be written as a linear combination of the system of polynomials (|3.2|) as follows 

n — m — 1 

£( Q <« [i?Q!J (z) = ^ 6 R , n , n -kX n -k L n - k (z) 

k— — m—l 
n-1 

^ flfl, n ,n-Jk [Q n -fc(«)] , (3.13) 



k—n—m 



Taking £^ Q,/3 ^[-] on the left hand side of the equality (|3.12p , we get 



Z^[RQ' n ]{z) = R{z)C^[Q l n {z)] + {l-z 2 )p{z)Q' n {z) 

= \ n R(z)L n (z) + (l-z 2 )p(z)Q' n (z). (3.14) 



From (fTT]) 



m+l 



R{z)L n (z) = Yl 

(z). (3.15) 



/e=— m— 1 
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Substituting (jXTU|) and (|3~T5|) in (|3~14"]) . we have 

m+l 



£( Q ^[i?Q;](z) = ^ (A» e fl ,„,„_ fc + d„,„_ fc ) L„_ fc (z). (3.16) 

k— — rn— 1 

As n > 2 (to + 1), we can assume that relation (|3.16[) is the expansion of the polynomial 
jC( a 'P> [R Q' n ] in terms of the polynomials L n . Now, identifying coefficients between (|3.13j) and 
(I3.16[) we have that 0R, n ,n-k = for all k = 1, . . . , n — m — 2, and we have the formulas (|3.3[l - 
(EU). □ 

4 Asymptotic behavior of the sequence {Q n }n>m and their 
zeros 

In this section we study the asymptotic behavior of the polynomials Q n and their zeros. Let us 
denote by || ■ \\ a ,p and || • || M the standard norms on the spaces L 2 (fx a ^) and L 2 (/j,) respectively. 
The following result is essential in the proof of the theorems in this and the next section. 

Theorem 4.1. Let n £ V m (a, P) where m g N. Then 

Qn{z) 



(a>j9) —r P m *(P,*), (4- 1 ) 

Pn (Z) n-HX 

uniformly on closed subsets of<C \ [—1, 1] where C = Cljjoo}. 

First, we state a preliminary lemma which follows from Theorems 26 and 29 of [6] §6.1]. 

Lemma 4.2. Let /i € P m (a,^) where m e N. 7/ {L n }nez + w i/ie sequence of monic orthogonal 
polynomials with respect to fi 

Ln(z) ~> 



f„ (z) n^oo 

uniformly on closed subsets ofC \ [—1, 1]. 
Proof o/ Theorem\Jl\ . From and (ED 



Q n (z)-L n (z) ^ f X„ \ P i-k\ z ) 



As A„ = — n(l + n + a + /?), for each fc fixed, k — 1, 2, . . . , m, 

lim = 1. (4.3) 

Let if be a closed subset of C \ [— 1 , 1] . From the interlacing property of the zeros of consecutive 
Jacobi polynomials on [—1,1], it easily follows that there exists a constant M* such that for all 

z e K 

'""<M fe <M*, fc = l,...,m, (4.4) 
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where M k — sup \z — x\ k and M* = max{Mi, M m }. 

x6[-l,l] 

From (|2.3p . it is not hard to see that there exist two monic polynomials of degree 4 (to — k) in 
the variable n, q[ a l^_ k ) (n) and g^'fm-fc) (")' sucn that 



P, 



(",/3) 



n—k 



a,0 



= 4 



fe_ m «l,4(m-fe) 



(n) 



o (Q ' /3) (n) 



Q,/3 



fc = 1, 2, . . . , m. 



Therefore, from the Cauchy-Bunyakovsky-Schwarz inequality and the extremal property of the 
monic orthogonal polynomials, for n sufficiently large, we get 



On n — k 



< 



< 



\L r 



\a,/3 



\fT n -k 



< 



Cl 



Ln\\^< 



F, 



a,f3 



T n -k 

Cl 2 m - k 



Cl 

T n -k 



1 ( 2,il~k)( n ) < Ci 2" I+1 



\ a (a ' P) (n 

\ ( ilMm-k)y 



(4.5) 



where 1 < k < to, ci = sup p(x) and r is the leading coefficient of p. Hence by (|4.2p . (|4.3p . 

x6[-l,l] 

(g3), and g3|) 



- L n (z) 



0, uniformly on closed subsets of C \ [—1, 1], 



and from Lemma 14721 the asymptotic formula (|4.1[) is established. 
Corollary 4.3. Let p e V m (a,(3) where to e N. Then 
1. 

lim Q n {z) 

n— >oo 

uniformly on compact subsets ofC \ [—1, 1] 



1 



□ 



(4.6) 



2. The set of accumulation points of the zeros of the sequence of polynomials {Q n }n>m *s [— 1, 1], 
i.e. 



p| Q {*:&(«) = ()} = [-1,1]. 



For eac/i n at least (n — to) zeros of Q n are contained on [—1,1]. 

Proof. The first part of the theorem is an immediate consequence of relations (|4.1[) and jTOj (8.21.9) 
and (4.21.6)]. 

To prove the second part, from (|3.9|) it easily follows that Q n has at least n — m zeros contained 
in [-1,1] (cf. (El §3.3]). 

The function in the right-hand side of (|4.1j) in Theorem 14. II is holomorphic and does not have 
zeros in C \ [—1, 1]. Let K be a closed subset of C \ [—1, 1], from the Rouche's theorem we have 
that for n large the polynomial Q n does not have zeros on K , i.e. the zeros of the sequence of 
polynomials {Q n }n>m can not accumulate outside [—1, 1]. 
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On the other hand, (|4.6p implies the weak star asymptotic of the zero counting measures of 

the polynomials {Q n }n>m (cf- [HI Ch. 2]). That is, if we associate to each Q n the measure 

^ \ 

fi n = — > 5 U . then d[i n (x) — > ^=^= (the equilibrium distribution on [—1,1]) in the 

Q n {uj)=0 

weak-* topology and this implies that the zeros of {Q n } n >m must be dense in [—1, 1]. □ 



5 Asymptotic behavior of the sequence {Q n }n>m and their 
zeros 

Some basic properties of the zeros of Q n follow from (jl.6p . For example, the multiplicity of the 
zeros of Q n is at most 3, a zero of multiplicity 3 is also a zero of L n and a zero of multiplicity 2 is 
a critical point of Q n . 

From the second part of Corollary 14.31 we get that Q' n has at least (n — m — 1) zeros of odd 
multiplicity on the open interval ] — 1, 1[. For m = 1 we have that 

Theorem 5.1. Under the same hypothesis of Theorem \ 3.1l if m = 1 the critical points of Q n 
interlace the zeros of L n . 

Proof. If rn = 1 then from (|3.9p the polynomial Q n has at least (to — 1) real zeros of odd order 
on ] — 1, 1[. But, Q n is a polynomial with real coefficients and degree n, consequently the zeros of 
Q n are real and simples. As Q n — Q n , from Rolle's theorem all the critical points of Q n are real, 
simple and (n — 2) of them are contained on ] — 1 , 1 [. 

Denote P{x) = (1 — x) a+1 {l + x)" +1 Q' n (x). As a and (3 are real numbers in general P is 
not a polynomial. Notice that P is a real-valued differentiable function on [—1,1]. Without loss 
of generality, suppose that there exists x G]l,oo[ such that P(x) = 0, as P(l) — from the 
Rolle's theorem there exists x' e]l,x[ such that P'(x') — 0. But, from (11.31) and (|1.6I) P'(x) = 
A n (l — x) a (1 + x)P L n (x) and all the critical points of P are contained in [-1,1]. Hence all the critical 
points of Q n are contained in ] — 1, 1[. Again using the Rolle's theorem, it is straightforward that 
the critical points of Q n interlace the zeros of L n - □ 

From Corollary 14.31 we have that the set of accumulation points of Q' n is [—1,1]. For m = 1, 
Theorem 15.11 gives that the critical points of Q n are in [—1,1], interlace the zeros of L n , and are 
simple. Numerical experiments also show this behavior for m > 1. We conjecture that this always 
is the case. 

For the proof of Theorems 11.11 and 1 1 . 2 1 we will use the following result. 
Lemma 5.2. Let /i G V m (a,fi), where m £ N and define for z £ C, A(z) = sup \z — x\ 

a;6[-l,l] 

and 8{z) = inf \z — x\.. If {Cn}n>m is a sequence of complex number with limit ( £ C and 
xe[-i,i] 

{Qn}n>m the sequence of monic orthogonal polynomials with respect to the pair (C^ 01 '^, fi) such 
that Qn(Cn) = 0, then: 

1. For every d > 1 there is a positive number Nd, such that {z £ C : Q n {z) = 0} C {z £ C : 
\z\ < A(C) + d} whenever n > Nd- 

2. //<5(C) > 2, the zeros of Q n can not accumulate on [—1,1] and for n sufficiently large they 
are simple. 

Proof. We already know that Q n (Cn) = and if Q n {z) = then Q n {z) — Q n {Cn)- From the 
Gauss-Lucas theorem (cf. [HI §2.1.3]), it is known that the critical points of Q n lie in the convex 
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hull of its zeros and from 2. of Corollary 14.31 the zeros of the polynomials {Q n }n>m accumulate on 
[— 1, 1]. Hence from the bisector theorem (see the proof of Theorcm l5.1l or [51 §5.5.7] ) \z\ < A(£„)+l 
and the first part of the theorem is established. 

To verify the second assertion of the theorem, note that if z is a zero of Q n , from ([2. 8(1 we get 



n 

k=i 



%7l,k 



(5.1) 



Let V e ([— 1, 1]) = {z £ C : 8{z) < e} be an e-neighborhood of [—1, 1]. On the other hand, as 
lim C„ = £, then for all e > there is a N e > such that |<5(Cn) — ^(C)l < £ whenever n> N e . 

n— >oo 

If 8(C) > 2, let us choose e = = - (5(Q — 2) and suppose that there is a zq £ V e( ([— 1, 1]) 
such that Qn(zo) = for some n > N e . Hence 



n 



%0 %n : k 



< 



2 + e c 

*(Cn) 



< 1, 



(5.2) 



which is in contradiction with (|5.1[) . Hence {z £ C : Q n (z) = 0} f] V £n ([— 1, 1]) = for all n > N e( , 
i.e. the zeros of Q n can not accumulate on V E( ([—1,1]). 

From ()2.8|) it is straightforward that a multiple zero of Q n is also a critical point of Q n - But, 
from 2. of Corollary 14.31 and the Gauss-Lucas theorem the critical point of Q n accumulate on 
[—1,1]. Thus, we have that for n sufficiently large the zeros of Q n are simple. □ 



Proof of Theorem From (|2.8|) the zeros of Q n satisfy the equation 



Qn(z) " = Qn(Cn) 



(5.3) 



If z £ C \ [—1,1], by taking limit when n oo, from 1. of lemma I5T21 and using (|4.6p on both 
sides of (|5.3p . we have that the zeros of the sequence of polynomials (Q n )n>m cannot accumulate 
outside the set 

{zeC: \z+Vz 2 -l\ =e"«] |J[-1,1]. 

Hence z + ^z 2 - 1 = e^ +lf) and z - V^ 2 - 1 = e'^i+M) for < 9 < 2tt, where we have that 
2z = e^ +ie + er^i+^y 

The assertion for S(() > 2 is straightforward from 2. of Lemma T5. 2 1 □ 
Now, we will state the relative asymptotic between the polynomials {Q n }n>m and the corre- 



sponding Jacobi polynomials P„ 



Proof of Theorem \1.2[ 
1.- Let us prove first that 



Qn(z) 



1 



1, 



(5.4) 



Qn(z) Qn(z) 

uniformly on compact subsets K of the set {z £ C : |<y?(z)| > |y(C)|}- I n order to prove (|5.4|) it is 
sufficient to show that 

Qn(Cn) 



Qn(z) 



0, uniformly on X. 



(5.5) 
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From [TQl (8.21.9) and (4.21.6)], we have the well known strong or power asymptotic of the monic 
Jacobi polynomials 



V n (z) ™ V 2(«-l)y V 2(« + l) 
uniformly on compact subsets of C \ [—1,1]. Note that 



(5.6) 



Qn(*) Pi^(Cn) Qn(z) ^ (Cn) 2"P^(z)WW 

From (|4.1[) and (|5.6[) the first four factors in the right hand side of the previous formula have finite 
limits; meanwhile, the last factor tends to when n — > oo, and we get (|5.5|) . Finally the assertion 
1 is straightforward from (14. lj) . 



2.- For the assertion 2 of the theorem it is sufficient to prove that 

Qn{z) Q n (z) 



1 



1. 



Qn(Cn) Qn (Cn) 

uniformly on compact subsets -KT of the set {z £ C : | | < |^(C)|} \ [ — 1, 1]- Note that 



(5.7) 



(z) P^\Cn)2 n P^\z) V n (Cn) f<p(z)Y 



Qn(C») " P^*) Qn(Cn) W 2™ P^ ((„) V ^(Cn) 

Now, the first part of the assertion 2 is straightforward from (|4.1I) . 

If (5(C) > 2, let V e ([-l,l]) = {z e C : (5(z) < e} be a e-neighborhood of [-1,1], where 

E-E - 5i0 



— 1. By the same reasoning that was deduced (15.21) we get that 



?«(*) 



,(C») 



< k™, for all z eV e ([-l, !]),«< 1. 



(5.8) 



Hence from the first part of the assertion 2 and 



we get the second part of the assertion 2. □ 



6 Fluid dynamics model of sources and stagnation points 

The fluid dynamic interpretation that we will consider in this section was introduced by H. Pijeira 
et al in [2] . In that paper the hydrodynamic model was a reinterpretation of the electrostatic model 
studied by H. Pijeira et al in [7]. The difference between the fluid dynamic model in [5] and the 
model introduced in the present paper is the complex potential used. 

Let us consider a flow of an incompressible fluid in the complex plane, due to a system of n — 1 
source points (n > 1) fixed at Wi, 1 < i < n — 1, with unitary rate of fluid emission per unit 
time (strength of the source), and two additional source points at 1 and —1 with strength a > 
and b > respectively. Here, a source is a point in which the fluid is continuously created and 
uniformly distributed in all directions with constant strength (steady source). Let us call flow field 
generated by a Jacobi set of sources to a flow of a fluid under the above conditions, or simple a 
flow field. 



12 



The complex potential of a flow field at any point z (cf. [31 Ch. 10] and [SJ Vol. Il-Ch. 6]), by 
the superposition principle of solutions, is given by 

n-l 

T(z) = Vlog(z- Wl )+alog(z-l) + 6 1og(z + l), 



»=1 



n-l 



log \{z-\) a (z + \) b Y[(z- Wi ) . (6.1) 



2 = 1 



From a complex potential T, a complex velocity V can be derived by differentiation (V(z) = 
^j(z)). A standard problem associated with the complex velocity is to find the zeros, that corre- 
spond to the set of stagnation points, i.e. points where the fluid has zero velocity. 

We are interested in an inverse problem in the following sense, build a flow field such that the 
stagnation points are at preassigned points with nice properties. As it is well known, the zeros of 
orthogonal polynomials with respect to a finite positive Borel measure on [— 1, 1] have a rich set of 
nice properties ( \10\ Chapter VI]), and will be taken as preassigned stagnation points. Here, we 
consider that fi <E Vi(a, (3). In the next paragraph the statement of the problem will be established. 

Problem. Let {x\, X2, ■ ■ • , x n } be the set of zeros of the nth orthogonal polynomial L n with 
respect to fi € Vi(a, f3) with 1 < n. Build a flow field (location of the source points w\, . . . , w n —\) 
such that the stagnation points are attained at the points X4, i = 1, 2, . . . , n. 

Let Q n be a monic polynomial of degree n, and denote its critical points by {w\, W2, ■ ■ ■ , w n _i}, 
thus 

n-l , v 

Q'M = nJliz-Wi), T (z) = log (- (z-l) a (z + if Q' n (z)), 

i=l V '' 



w , dT.. ((z-ir(z + l) b Q' n {z))' C^-V[Q n ]{z) 



dz w (z- 1)« (z + l) b Q' n (z) (z -l)(z + 1) Q' n (z) • 
dV 

From (16.11) and Theorem 15.11 -7— (a;^) = for each stagnation point (zeros of L n ), k = 

dz 

l,2,...,n, i.e. 

£ (o - 1 ' 6 - 1) [Q„](x fe ) = 0, fc = l,2,---,n. (6.2) 



From Corollary 12.21 there exists a monic polynomial Q n of degree n, unique up to an additive 
constant, satisfying equation (|6.2I) . i.e. 

C {a - 1 > b - 1 \Q n ]{z) = \ n L n {z) 1 X n = -n(n+a + b-l). (6.3) 

Note that (|6.3p is the same as (|1.6p with a = a — 1 and j3 = b — 1. Therefore the n — 1 source 
point of the flow field {u>i, . . . , w n —i} are the critical point of the nth orthogonal polynomial with 
respect to the differential operator £( a_1 ' fc_1 ). 

Answer. A flow fields generated by a Jacobi set of sources with complex potential (|6.ip and 
preassigned stagnation points at the zeros of the nth orthogonal polynomial with respect to the 
measure fi € Vi(a,/3) with n > 1, has its sources points (with unitary strength) located at the 
critical points of the nth orthogonal polynomial with respect to (£S a '^\ /i). 

In Theorem l5.11 we proved that for m = 1 all the critical points of Q n are simple, contained in 
[—1,1] and interlace the zeros of L n . At the beginning of the Section [SJ we conjectured that this 
theorem is true for all m £ N. If this were true, then it is not difficult to see that the above model 
holds for m 6 N. 

Note that, if we consider a system of electrostatic charges with potential given by (|6.1|) instead of 
a system of source points with the same potential function, then we have an analogous electrostatic 
interpretation. 
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As it is known, the zeros of the Jacobi polynomials have an electrostatic interpretation (see 
[TU1 §6.7]) as the equilibrium points of a certain potential function. For the case of orthogonality 
with respect to a differential operator the electrostatic interpretation is an inverse problem in the 
sense that the equilibrium points are known and the question is to build the electrostatic field. 
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